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■ Abstract 

(N 

We present an empirical method that uses multicolor light curve shapes (MLCS) 

> 

' to estimate the luminosity, distance, and total line-of-sight extinction of Type la su- 

pernovae (SN la). The empirical correlation between the MLCS and the luminosity 
, is derived from a "training set" of nine SN la light curves with independent distance 

' and reddening estimates. We find that intrinsically dim SN la are redder and have 

Q^' faster light curves than the bright ones which are slow and blue. By thirty-five days 



after maximum the intrinsic color variations become negligable. A formal treatment 



^ \ of extinction employing Bayes' theorem is used to estimate the best value and its un- 



certainty. Applying MLCS to both light curves and to color curves provides enough 
information to determine which supernovae are dim because they are distant, which 
are intrinsically dim, and which are dim because of extinction by dust. The precision 
of the MLCS distances is examined by constructing a Hubble diagram with an inde- 
pendent set of twenty SN la's. The dispersion of 0.12 mag indicates a typical distance 
accuracy of 5% for a single object, and the intercept yields a Hubble constant on the 
Cepheid distance scale (Sandage et al 1994, 1996) of Hq = 65± 3 (statistical) km s~^ 
Mpc^^ (±6 total error). The slope of 0.2010 ± 0.0035 mag over the distance interval 
32.2 < fi < 38.3 yields the most precise confirmation of the linearity of the Hubble law. 

subject headings: supernovae: general ; cosmology: distance scale 
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1 Inhomogeneity of Type la Supernovae 

Since the Curtis- Shapley debate of 1920 (Curtis 1921, Shapley 1921), the determination 
of supernova (SN) luminosities has been central to the discussion of cxtragalactic distances. 
Shapley (1919) argued against the "Island Universe" hypothesis, because it required certain 
novae (such as S Andromedae of 1885) to reach the astonishing luminosity of M = —16. He 
considered this to be "out of the question". Curtis (1921) countered, concluding that "the 
dispersion of the novae in spirals and in our galaxy may reach ten absolute magnitudes... a 
division into two classes is not impossible" . This distinction between novae and supernovae, 
required by the extragalactic nature of the nebulae, was later made explicit by Baade and 
Zwicky (1934). They showed that in addition to their tremendous difference in absolute 
luminosity, the photometric and spectroscopic behavior of supernovae is distinct from novae. 
Baade (1938) showed that supernovae were more uniform than novae, with a dispersion at 
peak of 1.1 magnitudes, making them suitable as extragalactic distance indicators. 

The precision of supernova distance estimates has increased as the SN la class has been 
better understood and more narrowly defined. The low dispersion in Baade's sample bene- 
fited from the fortuitous absence of Type II supernovae, which are significantly less luminous 
in the mean. Beginning with SN 1940B, Type II supernovae were classified by the presence 
of hydrogen in their spectra (Minkowski 1941). The growing hst of spectroscopically defined 
Type 1 supernovae had dispersions at peak of 0.8-0.6 magnitudes (Minkowski 1964, Kowal 
1968, Kirshner et al. 1973, Oke and Searle 1974). However, this sample included a number 
of "pecuhar" SN la noted for their lack of sihcon, which are now recognized to arise from 
massive stars that lose their envelope before core collapse (see Wheeler & Harkness 1990). 
After removing these silicon-deficient objects, now classified as Ib's and Ic's (Doggett and 
Branch 1985, Uomoto and Kirshner 1985, Wheeler and Levreault 1985, Wheeler and Hark- 
ness 1986), the remaining Type la supernovae (SN la) form a more homogeneous set which 
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serve as even more precise indicators of astronomical distances. Leibundgut (1989) devised a 
set of standard templates to describe the photometric behavior of SN la and to estimate the 
peak apparent magnitude. Hubble diagrams constructed using the peak of photographic SN 
la light curves had observed dispersions ranging from aM— 0.65-0.36 magnitudes depending 
on which objects and color bands were used (Tammann and Leibundgut 1990, Branch and 
Miller 1993, Miller and Branch 1990, Delia Valle and Panagia 1992, Rood 1994, Sandage and 
Tammann 1993, Sandage et al 1992,1993,1994). An ambitious program to cahbrate nearby 
SN la through Cepheid variables observed with the Hubble Space Telescope has been under- 
taken by Sandage et al (1992, 1993, 1994, 1996). Assuming SN la to be homogeneous yields 
a Hubble constant in the range 50-58 km s"^ Mpc"^ (Sandage et al 1992,1993,1994, 1995, 
Schaefer 1994, 1995a,b, 1996, Branch & Tammann 1992) with the most recent measurement 
giving 57 ± 4 km s~^ Mpc~^. We show in §6 and §7 that the precision of this measurement 
is improved and the value of Hq altered by including information contained in the light and 
color curve shapes. 

The hypothesis that SN la are standard candles drew support not only from empirical 
studies, but also from the earliest theoretical models which suggested they arose from igni- 
tion of a carbon-oxygen white dwarf at the Chandrasekhar mass (Hoyle and Fowler 1960, 
Arnett 1969, Colgate & McKee 1969). In these models a supersonic shock wave travels 
through the degenerate star, burning material into ^^Ni at a temperature of five billion 
degrees (Khokhlov, MuUer, and Hoflich 1993, Mazurek & Wheeler 1980). Because the det- 
onation is supersonic, the pre-shock region cannot expand to decrease the pressure or the 
burning temperature. Further, the Fermi pressure of the degenerate material in the post- 
shock region remains insensitive to temperature for longer than the burning timescale. The 
result is a total incineration and the production of a pure mass of nickel. Such a standard 
explosion of a uniform mass would lead to a homogeneous light curve and uniform luminosity. 
Yet, these complete burning models of la's do not reproduce the intermediate mass elements 
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which are seen in the spectra of SN la (Wheeler and Harkness, 1990). A successful model 
(Nomoto, Thielemann, and Yokoi 1984) which matched the observational constraints was so 
persuasive that Arnett, Branch, and Wheeler (1985) and Branch (1992) suggested calibra- 
tion of the Hubble constant based only on theoretical models of uniform nickel production. 
However, a variety of models (Livne 1990, Khokhlov, MuUer, and Hofiich 1993, Woosley 
& Weaver 1992, Hofiich, Khokhlov, & Wheeler 1995) match the observed features of the 
spectra and produce a range of nickel masses, and a range of predicted luminosities. These 
models employ subsonic deflagration fronts, pulsations, or off-center explosions to allow the 
surface layers to pre-expand and burn at low temperature. The success of these models in 
reproducing the observed spectra opens a large range of theoretical possibilities. Unlike the 
first monoenergetic models, these models suggest a wide range of luminosities might result 
from the ignition of a white dwarf. 

Recently, precise observations of SN la made with CCD detectors show evidence for 
inhomogeneity in both luminosity and light curve shape. One of the first SN la observed 
with a distinctly different light curve was 1986G (Phillips ct al 1987, Cristiani et al 1992) 
which displayed a spectacularly rapid decline in its B and V light curve and unique spectral 
characteristics including stronger-than-usual Si features. Although SN 1986G was heavily 
reddened by dust, reddening cannot significantly alter the shape of the light curve. SN 
1991bg in NGC 4374 is the most extreme SN la in an increasingly apparent photometric and 
spectral sequence. Leibundgut et al (1993) (also Fihppenko et al 1992) described a number of 
photometric abnormalities of 1991bg with respect to his templates. These include the fastest 
post-maximum decline (2.05 and 1.42 mag decrease drop in B and V in the fifteeen days after 
maximum compared to 1.22 and 0.64 mag for the templates in B and V), a narrow luminosity 
peak, and an intrinsic red color near maximum. A simple and convincing argument that SN 
la have a large spread in luminosity is that SN 1957B, which occurred in the same galaxy, 
was 2.5 magnitudes brighter in B than SN 1991bg. In addition, SN 1991bg was at least 2 
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magnitudes fainter than other SN la in the Virgo cluster, of which NGC 4374 is a member. 
This extreme SN la seems to have a twin in the sub-luminous SN 1992K (Hamuy et al 1994), 
which strongly resembles 1991bg in photometric and spectral behavior. 

At the opposite extreme of the SN la class, SN 1991T showed spectral and photometric 
peculiarities which were different from those seen in the rapid decliners. Phillips et al (1992) 
found the light curves in B and V to rise and decline more slowly than the standard templates 
near maximum, and a month after peak, this shallower decline resulted in a light curve 0.2-0.3 
mag brighter than the templates. Although SN 1991T's host galaxy, NGC 4527, appears to 
lie south of the main Virgo cluster, Phillips et al (1992, 1993) estimates the peak luminosity 
exceeded that of other SN la in Virgo by 0.3-0.5 mag. Prom the narrow Na I D absorption 
hne, Fihppenko et al (1992) deduced that SN 1991T is dimmed by dust in NGC 4527 and 
concluded it may have been as much as ~ 0.9 mag more luminous than a typical SN la. 
More recently, the Calan/CTIO supernova search yielded SN 1992bc and 1992bo (Maza et 
al 1994), two SN la with similar recession velocities of 6050 and 5660 km s"^, but with peak 
apparent luminosities differing by 0.69 mag in B. The large difference in apparent magnitude 
and the small difference in recession velocity imply that SN 1992bc was intrinsically brighter 
by 0.8 ± 0.2 mag than 1992bo. Por a difference in distance to account for this difference in 
luminosity, the peculiar velocities of the two supernovae would have to differ by ~ 2,500 
km sec~^, an unlikely alternative. Interestingly, SN 1992bc declined more slowly than the 
average template while SN 1992bo's post-maximum fall was more rapid than the average 
template, a result in good accord with later analysis by Phillips (1993). 

Even before these recent examples of inhomogeneity, less precise photographic measure- 
ments by Barbon et al (1973) suggested that there might exist two photometric classes; those 
with "fast" decline rates after maximum which were intrinsically brighter supernovae, and 
those with "slow" decline rates which were fainter. With the poor quality of photographic 
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and photoelectric photometry then available, such a real distinction was difficult to demon- 
strate convincingly. Pskovskii (1977,1984) suggested a continuous photometric sequence of 
SN la light curves. He introduced a parameter, (3: the slope of the B band post- maximum 
decline in magnitudes per 100 days. Using 54 photographic SN la light curves, Psovskii 
found a weak correlation between /? and the absolute B magnitude at maximum light which 
was opposite to that of Barbon (1973). Early difficulties in measuring a relation between SN 
la light curve shape and intrinsic luminosity resulted from noisy photographic data which 
were poorly sampled and calibrated. These difficulties were compounded by the problem of 
measuring a decaying light curve on a bright galaxy with a non-linear photographic detector 
(Boisseau and Wheeler 1991). With the advantage of better data measured with linear de- 
tectors, Phillips (1993) demonstrated conclusive evidence for a luminosity-light curve decline 
relation among SN la. Using a set of well-sampled SN la hght curves with precise optical 
photometry and accurate relative distances, Phillips found that the absolute luminosity in 
B,V, and I is correlated with the B band decline in the fifteen days following maximum light. 
The sense of the correlation is that dimmer SN la fall more rapidly after maximum than 
bright SN la. Application of this relation to his sample results in a significant reduction of 
the dispersion in B,V, and I luminosity from 0.8, 0.6, and 0.5 mag to 0.36, 0.28, and 0.38 re- 
spectively. A more recent investigation by Tammann and Sandage (1995) has examined the 
luminosity-light curve decline relation among "normal" SN la, in this case "normal" is de- 
fined as having a {B — V)max < 0.30 mag. They tentatively confirm a luminosity dependence 
on light curve decline, with a slope that is shallower than Phillips (1993) but consistent with 
Hamuy et al (1995) for a similar set of "normal" SN la. In §6 we show for a sample restricted 
to "normal" SN la, a significant correlation between light curve shape and luminosity exists. 

This empirical relation between light curves and luminosity in SN la has been paced by 
an abundance of theoretical models which can account for the observed behavior (Hofiich, 
Khokhlov, and MuUer 1993, Woosley and Weaver 1994, Ruiz-Lapuente et al 1993, Livne 
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and Arnett 1995, Hoflich and Khokhlov 1995). These new models include deflagration burn- 
ing fronts, off-center detonations, surface helium burning, pulsed delayed detonations and 
sub-Chandrasekhar progenitors. These models give plausible causes for the observed inho- 
mogeneity of SN la and for the origin of the empirical relations between light curve shape 
and luminosity. 

Most recently, Hamuy et al (1995) have employed a template- fitting approach and we 
(Riess, Press, and Kirshner 1995a, hereafter RPK 95a) have developed a hnear estimation 
algorithm to use the distance-independent light curve shapes to improve the precision of 
distance measurements to SN la. The techniques have much in common and both yield 
Hubble diagram dispersions of ~ 0.2 magnitudes for an overlapping set of objects. The 
Light Curve Shape (LCS) method, described in RPK 95a, has the advantage of providing 
quantitative error estimates for the distance. The present paper extends the LCS technique 
to use the shapes of B— V, V— R, and V— I color curves which provide enough information 
to determine the relation between absolute luminosity and intrinsic color. Knowledge of 
an SN la's intrinsic color allows us to measure the extinction from the observed reddening. 
For each well-observed SN la we measure the luminosity, extinction and the extinction- 
corrected distance. The multicolor light curve shape (MLCS) method significantly increases 
the precision of distance estimates from SN la hght curves as we show in §6. 

There are many potential applications for a bright distance indicator with < 10% pre- 
cision. Nearby (0.01 < z < 0.1), it should be possible to measure the Hubble constant to an 
accuracy which is limited only by the underlying calibration of Cepheids. It is important to 
compare the Hubble constant derived using the light curve shape-luminosity relation with 
determinations which have assumed a homogeneous luminosity for SN la (Sandage et al. 
1992, 1993, 1994, 1996). Using the velocity residuals from Hubble flow, we have measured 
the motion of the Local Group with respect to the rest frame of galaxies with supernovae 
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(Riess, Press, and Kirshner 1995b). At even greater distances (0.3 < z < 0.6) MLCS dis- 
tance measurements of all well-observed SN la could be used to determine the cosmological 
deceleration parameter, Qo (Perlmutter et al 1995, Perlmutter et al 1996, Schmidt et al 1996, 
lAUC 6160). 

Some have sought to improve the homogeneity of the observed SN la by restricting the 
sample to supernovae with "normal" spectra (Branch, Fisher, and Nugent 1993, Sandage 
et al 1994, 1996). While spectroscopic information may prove useful in producing smaller 
dispersion in distance estimates, the difficulty in obtaining good spectra of very distant SN 
la makes it hard to identify subtle spectral variations. A sample cut which cannot be applied 
with equal effectiveness for nearby and for distant SN la could lead to a bias in cosmological 
parameters determined by them. We prefer to develop a method that can be applied to all 
SN la. 

Given at least one light curve and one color curve, photoelectrically observed within 
ten days after maximum, our MLCS method can distinguish between the effects of distance, 
intrinsic dimness, and dust for all SN la. In §2-5 we develop the multicolor hght curve shape 
method for measuring extinction-corrected distances. In §6, for an independent sample, 
we compare the extinction-corrected MLCS distances with distances determined by more 
limited assumptions. In §7, we estimate the Hubble constant and discuss sample membership, 
selection bias and the range of progenitors which may be responsible for the empirical range 
of SN la luminosities and colors 

2 Leeirning Curves 

SN la light curves are noisy unevenly sampled time series, which, when diligently ob- 
served, are available in four bands, B, V, R, and I. Although some SN la are found in 
ellipticals, most are found in spirals and irregulars where they can be subject to significant 
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extinction by dust. At present, the number of well sampled SN la light curves available on a 
modern photometric system is limited (i.e. < 50) (Hamuy et al 1995, Riess et al 1996, Ford 
et al 1993). 

We employ a general model to establish the empirical relation between light and color 
curve shape and luminosity from a nearby subset of SN la with accurately known relative 
distances and extinctions. Then we use these results to estimate the extinction-corrected 
distance for a separate distant sample solely from the observed light and color curves. The 
observables are a visual band light curve and up to three color curves. 

To first approximation, the observed light and color curves of SN la are homogeneous, 
resembling standard template light and color curves with the addition of noise (n) and 
individual offsets that result from the apparent distance modulus (//) and a color excess 
{E color) due to reddening by dust: 

mv = Mv + A* + n (1) 
niB - mv = (B - V)o + Eb-v + n (2) 
mv - niR = (V - R)o + Ey-R + n (3) 

mv - mi = (V - I)o + Ev-i + n (4) 

Here bold-faced denotes a vector whose entries are measured or determined as a function 
of time. In standard convention, M is an absolute magnitude, m is an apparent magnitude, 
and {B — V)q, {V — R)q, and {V — J)o are unreddened color curves. Each available color 
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curve yields an Ecoior term which, combined with a standard extinction curve, reduces the 
uncertainty in measuring the total line-of-sight extinction. Ay- This description of light and 
color curves assumes intrinsic homogeneity of SN la, and is equivalent to the "template" 
approach of Leibundgut (1989), Sandage et al (1994, 1996), and Sandage and Tammann 
(1993). 

Improvements in the quality of the available data set and the success of Phillips (1993) 
motivate an approach which, to first order, can account for the observed variations in the 
light and color curves and intrinsic luminosity. The most economical approach is to adopt 
a single parameter and correlate it with the variations of the observed curves. The natural 
parameter to choose is the amount by which the intrinsic luminosity differs from an SN la 
of standard brightness and curve shape since, in the end, it is this difference we hope to 
measure. We call this the "luminosity correction". A, where A = My — My^standard, and 
Mv,standard IS the lumiuosity of the SN la chosen to depict the "standard" SN la event. This 
difference is measured by convention at the date of B maximum. The functions of time 
which, to first order in A, correct the observed curves to the template curves are "correction 
templates", Rv{t) or Rcoior{t)- Using these definitions, we expect two supernovae whose 
absolute visual luminosities differ by A magnitudes to have light and color curves which 
differ by Rv{t)A and Rcoior{t)A. The assumption that luminosity correlates linearly with 
light and color curve shape need not provide a complete description, but given the limitations 
of our real data sets, it is a reasonable way to begin. The evidence that this is a useful way 
to proceed is provided in §6 where we show that this method produces a significant decrease 
in the dispersion for an independent set of SN la in the Hubble fiow. 

Our improved model for the light and color curves is thus 

mv = Mv + A* + RvA + n (5) 

where Rv{t)A. gives the deviation from the template, Mv{t), at each time. Similarly, 
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me - mv = (B - V)q + Eb-v + Rb-vA + n 



(6) 



mv - niR = (V - R)o + Ey-R + Rv-rA + n (7) 

mv - mi = (V - I)o + Ev-i + Rv-iA + n (8) 

where {B — V)o, {V — R)o, and {V — I)o are the unreddened template color curves. 
It is instructive to compare equations (5) through (8) which allow SN la inhomogeneity to 
equations (1) through (4) which impose SN la homogeneity. 

In matrix format we can write equations (5) through (8) as one system of equations 
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(9) 



where we have combined the Emior terms into one simultaneous measurement of Ay with the 



aid of the standard reddening law parameterization, ^^^^^^ = 3.1, ^q/Lr) — 3-9, E{\r-i) ~ 
1.9 (Savage & Mathis, 1979). Here, we assume the Galactic reddening ratios are vahd for the 
dust in distant galaxies. Elsewhere, we examine the homogeneity of these ratios in distant 
galaxies, and find the most likely form of the reddening law for dust in distant galaxies is 
consistent with the Galactic reddening law (Riess, Press, Kirshner 1996b). 
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The of the fit between data and model is 

= (n)^C-^(n) (10) 

where n is the vector of residuals in equations (5) through (8) or equivalently in equation 
(9). Here, C is the correlation matrix whose elements are intended to reflect the errors 
in observations and residual, unmodeled correlations in SN la light and color curves. The 
correlation matrix is discussed in detail in §4. 

Rybicki and Press (1992) have derived two analytical minimizations of the given in 
equation (10). One gives the best estimate of the correction templates, Ryit) or Rcoior{t), 
provided the SN la parameters, [/x Ay A] , are known, the other gives the best estimate of the 
parameters provided the correction templates are known. We employ both, using the nearby 
set with accurate relative distances as a "training set" to estimate Rv{t) and Rcoiorif)- Once 
trained, we measure the yu, Av, and A parameters for an independent set of SN la from the 
shapes of their light and color curves. 

Analytical minimization of in equation (10) with respect to Rv{t) or Rcoiorif) is 
independent of C and gives; 



or 



rCcolor[t) - ^^^j (i^j 

where the angle brackets denote an average over the training set weighted by the uncertainty 
in the given estimates of A. 

It is convenient to choose a particular object as the "standard" SN la whose light and 
color curves deflne the standard template curves, Mv{t) and Mcoiorif)- By definition, the 
"standard" SN la have a peak luminosity variation of A =0. The utihty of the method 
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Table 1: Training Set 



SN 


Galaxy 




Method 


IxRef 


E{B-V){a) 


Mv 


A(a) 


Phot.Ref 


1980N 


N1316 


31.15(0.10) 


SBF,PNLF 


2 


0.00(0.02) 


-18.71 


-0.17(0.10) 


1 


1981B 


N4536 


30.50(0.30) 


T-F 


1 


0.00(0.02) 


-18.54 


0.00(0.32) 


2 


1986G 


N5128 


27.72(0.10) 


SBF,PNLF 


2 


0.60(0.10) 


-18.13 


0.41(0.31) 


3 


1989B 


N3627 


29.40(0.30) 


T-F 


1 


0.35(0.03) 


-18.50 


0.04(0.31) 


4 


1990N 


N4639 


31.40(0.30) 


T-F 


1 


0.01(0.02) 


-18.82 


-0.28(0.31) 


5 


1991T 


N4527 


30.60(0.30) 


T-F 


1 


0.00(0.02) 


-19.10 


-0.56(0.31) 


6,7 


1991bg 


N4374 


31.05(0.10) 


SBF,PNLF 


2 


0.00(0.02) 


-17.10 


1.44(0.10) 


8,9 


1992A 


N1380 


30.65(0.11) 


SBF 


4 


0.00(0.02) 


-18.10 


0.44(0.13) 


10 


1994ae 


N3370 


31.28(0.30) 


T-F 


5 


0.14(0.03) 


-18.65 


-0.13(0.30) 


11 



11 references- (1) Pierce 1994; (2) CiarduUo, Jacoby, &; Tonry 1993; 
(4) Tonry 1991 (5) Dell Antonio 1995 

photometry references- (1) Hamuy et al 1991;(2) Buta and Turner 1983; (3) Phillips et al 1987 
(4) Wells et al 1993; (5) Leibundgut et al 1991; (6) Phillips et al 1992; (7) Ford et al 1993; 
(8) Filippenko et al 1992;(9) Leibundgut et al 1993; (10) Suntzeff et al 1996 (11) Riess et al 1996 



is insensitive to the choice of the template curves and luminosity since all quantities are 
defined relative to it. Leibundgut's templates, made from SN 1989B, SN 1980N, and SN 
1981B approximate the fight curves of the "normal" training set supernovae in B and V. For 
R and 1, we have constructed our own templates from SN 1989B, SN 1980N, and SN 1981B. 

In solving for the correction templates, we only need consistent relative distances for 
our training set of SN la since we correlate luminosity differences with light curve shape 
variation. No choice of absolute distance scale is necessary since all luminosity corrections, 
A, are relative to the standard SN la. Surface Brightness Fluctuations (SBF), Planetary 
Nebula Luminosity Function (PNLF), and TuUy and Fisher's luminosity-line width relation 
(T-F) provide accurate and consistent bias-corrected relative distances to the host galaxies of 
the training set of SN la (Strauss & Wilfick 1995, Jacoby et al 1992, CiarduUo, Jacoby, and 
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Tonry 1993, Pierce 1994, Kennicutt, Preedman, and Mould 1995). We have no dogma about 
the correctness of this absolute distance scale-in §7 we calibrate the absolute luminosity of 
a standard SN la, contained in Mv{t) with the Cepheid observations from Sandage et al 
(1994, 1996). This approach has the advantage of placing our distance scale directly onto 
that of the Cepheid variables. 

The training set of supernovae we have used also has precise optical photometry, well- 
sampled hght curves, and estimates of E(B — V) (see table 1). The color excesses in table 
1 arc estimated from comparisons with unreddened, photometrically similar SN la (Phillips 
1993, Wells et al 1994) and are used to correct the luminosity and color curves of the affected 
objects. In §7 we discuss our training set membership. On the SBP-PNLP-TP distance 
scale, SN la which define the template have an < My >=-18.54 on the date of B maximum; 
therefore each supernova's A = My — (—18.54). We have calculated the correction templates 
from equations (11) and (12) using the training set of supernovae provided in table 1. The 
correction templates and standard templates are listed in table 2. The standard template, 
Mv{t), as listed in table 2, is calibrated on the Cepheid distance scale discussed in §7. 

Adding various amounts of the correction templates to the standard templates generates 
an empirical family of light and color curves (see figure 1 ) . This family of curves demonstrates 
some interesting relations between light curve behavior and luminosity. The natural history 
of SN la is that intrinsically dim SN la rise and fall rapidly in V as compared to the leisurely 
rise and fall of intrinsically bright SN la. These results are similar to those of Phillips (1993) 
with the advantage that they show the SN la behavior from before maximum and more than 
fifteen days after maximum. 
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Table 2: Training Vectors 



phase standard templates correction templates weighting vectors 





Mv(t) 




Mv-R(t) 


My_,(t) 


Rv(t) 


RB-v(t) 


Rv-rW 


Rv-i(t) 








,Tj(t) 


-11.000 


-18.061 


-0.307 


0.007 


0.390 


1.257 


0.492 


0.273 


0.408 


0.000 


0.371 


0.130 


0.251 


-10.000 


-18.374 


-0.326 


0.009 


0.294 


1.257 


0.493 


0.273 


0.408 


0.002 


0.366 


0.112 


0.221 


-9.000 


-18.534 


-0.285 


0.012 


0.260 


1.267 


0.493 


0.273 


0.408 


0.008 


0.349 


0.112 
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Figure 1: Empirical family of SN la light and color curves parametrized by luminosity. This 
sample of empirical V, B-V, V-R, and V-I curves is derived from the training set and depicts 
the entire range of light and color curve shapes and their correlation with luminosity (on the 
Ccpheid distance scale). This set is obtained by adding the correction templates, i?y(t) or 
Rcoior{t), multiplied by various luminosity corrections, A, to the standard templates to make 
the best reconstruction of an SN la light and color curves. Intrinsically dim SN la rise and 
fall faster in V and have redder colors before day 35 than intrinsically bright SN la. After 
day 35, all SN la have more uniform colors. Prom the multicolor hght curve shape (MLCS) 
method we estimate the luminosity and extinction by dust independently from the distance 
to measure the extinction- free distance. Data shown as reconstructed, 91T=o, 94ae=n, 
86G=X,91bg=o, 92A=+, 80N=A 
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An interesting new result is the quantitative relation between SN la luminosity and 
color. Supernova luminosities correlate with their intrinsic colors at early times (see also 
Lira, P. 1995). Before day thirty-five, the dim SN la are red and the bright ones are blue. 
For example, at maximum, a decrease of 0.10 magnitudes in visual absolute luminosity 
corresponds to an increase in color (toward the red) of 0.05, 0.03, and 0.04 magnitudes 
in B— V, V— R, and V— I, respectively. Presuming that all SN la have a uniform color 
at maximum (as has often been done) is a poor assumption that can lead to incorrect 
predictions as we show elsewhere (Riess, Press, Kirshner 1996b). In particular, the relation 
between intrinsic color and luminosity helps explain how absorption "corrections" have led 
to increased dispersion in Hubble diagrams and to unphysical properties derived for dust in 
distant galaxies (Branch & Tammann 1992, Capaccioh et al 1990, Joeveer 1982, Tammann 
1987, Sandage et al 1993). For comparison, the color variations at maximum due to reddening 
for 0.10 magnitudes of visual extinction are expected to be 0.03,0.03, and 0.05 in B— V, 
V— R, and V— I. Although the sense of the color changes from absorption is the same as from 
intrinsic variation, the values are not. 

van den Bergh (1995) noted that the relation between intrinsic B-V color and lumi- 
nosity predicted by the mean behavior of Hoflich and Khokhlov's (1996) theoretical models 
coincidentally agrees with the standard reddening law. He takes advantage of this agreement 
to define a reddening-free luminosity which employs a single measurement of B-V color to 
account for both the luminosity variation intrinsic to the supernova and that which results 
from absorption by dust. Our empirical relation between luminosity and color shows that 
dust and intrinsic luminosity variation do not cause exactly the same change in color. This 
difference may cause the increase in dispersion around the Hubble line observed after van 
den Bergh's prescription is apphed (Riess, Press, Kirshner 1996b). 

More than thirty-five days past maximum light, all supernovae exhibit nearly uniform 
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colors. Entries in the correction templates for day thirty-five give differences of 0.00, 0.01, 
and 0.01 magnitudes in B— V, V— R, and V— I color for a 0.10 magnitude change in visual 
absolute luminosity. Detailed theoretical modeling of SN la show a similar relation between 
supernova color and luminosity (Hoflich & Khokhlov 1996). 

An alternate view of the photometric differences between intrinsically bright and dim SN 
la is presented in figure 2. The family of absolute B,V,R, and I fight curves show recognizable 
morphological variations. The B light curve family is similar in behavior to the V family; 
dim SN la rise and fall more rapidly in B and V than bright SN la. In R the brighter SN 
la have a "shoulder" ~ 25 days after B maximum. For dimmer SN la, this shoulder is less 
pronounced and disappears completely for the most underluminous objects. In the I band, 
the brighter SN la have two maxima. The first occurs quite early, ~ 5 days before the B 
band maximum. The second, broad maximum is at ~ 30 days after B maximum. As the 
luminosity of the SN la decreases, a number of changes in the I band light curve are apparent: 
the first maximum is later and broader while the second maximum is dimmer and occurs 
earlier. For the most underluminous SN la, the two maxima merge into one maximum which 
is broad and occurs ~ 5 days after B maximum. 
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Figure 2: Empirical family of SN la BVRI light curves parameterized by luminosity. This 
family of light curves is derived in the same way as the families in figure 1 and shows the 
differences in photometric behavior for bright and dim SN la. Intrinsically dim SN la rise 
and fall faster in B and V than intrinsically bright SN la. For the R light curve, a "shoulder" 
occurs ~ 25 days after B maximum in the bright SN la. This shoulder is weaker for dimmer 
SN la and is absent for the most underluminous ones. In the I band, the bright SN la 
have two maxima; one early (~ 5 days before B maximum) and one later (~ 30 days after 
B maximum). As the luminosity of the SN la decreases the first maximum occurs later 
and is broader while the second maximum is dimmer and occurs earlier. For the most 
underluminous SN la, the two maxima merge into one maximum which is broad and occurs 
~ 5 days after B maximum. Data shown as reconstructed, 91T=o, 94ae=n, 86G=X,91bg=o, 
92A=+, 80N=A 
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The light and color curve reconstructions in figure 1 provide a powerful means to measure 
the extinction-corrected distance at every phase of supernova observation. SN la appear dim 
because they are distant, obscured by dust, or intrinsically dim. We can distinguish between 
these possibilities by using the light and color curve shapes (which determine the intrinsic 
luminosity and color) and measuring the observed offsets (which determine the extinction- 
corrected distance). 

3 Extinction-corrected Distances from Multicolor Light Curve Shapes 

Figure 1 provides a ready guide for measuring extinction-corrected distances. Given 
BVRI light curve photometry, we seek the best set of curves for a fixed value of A which 
minimizes the between model and data. The offsets between the model and the data 
provide the best estimates of and Ay- Rather than searching the parameter space for 
a solution we take advantage of the linearity of our model for an immediate solution. 

Referring to the matrix form of the model in equation (9), we use the following defi- 
nitions; y is the column of apparent magnitude measurements, s is the column of standard 

templates, L is the three column matrix with correction templates and offsets, and q is the 
three element column of parameters. With these definitions, we rewrite equation (10) 
as 

= (y - s - Lq)^C-^(y - s - Lq). (13) 
The analytical minimization of with respect to the column of free parameters q gives 

q,,,, = Ay A]L, = (L^C-iL)-^L^C-^[y - s] (14) 

Equation (14) simultaneously racasmes the distance and the extinction using all the available 
light curve observations, reducing the reliance on any particular time of the supernova light 
curve. According to our definitions in equation (1) or (5), //is the apparent distance modulus 
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uncorrected for extinction. The extinction-corrected MLCS distance is given by ij, — Av. The 
"standard candle" distance, i.e. without any correction for hght curve shape (luminosity) 
or reddening, is given by /i + A. This is the distance derived by correctly fitting the shape 
of the visual light curve to find the peak, ignoring the color excesses, and assigning the 
SN la a standard luminosity. The standard errors for the parameters in q^est are given by 
the covariance matrix, (L"^C~^L)~^. These errors are the fitting errors which reflect the 
uncertainty in locating a light curve's best placement in figure 1, due to the presence of noise 
in either the training set light curves or in the independent light curve we are trying to fit. 

This approach, until now, treats the correction templates derived from our training set 
as if they were perfect. This is certainly not the case since the independent distance and Ay 
estimates for the objects that make up the training set objects are themselves not perfect. 
The same type of uncertainty is more easily seen in the Phillips (1993) relation where the 
light curve decline in the first fifteen days after maximum is correlated against the luminosity 
derived from independent distance and Ay estimates. The uncertainty in the slope of this 
relation must be considered when it is used. Fortunately, our correction templates are well 
constrained by the accuracy of the distance and Ay estimates combined with the size of 
the training set. The external source of error on the parameters can be found by varying 
the training set distances and Ay estimates in a Monte Carlo simulation to determine the 
effect on the fitted parameters. As expected, the external distance error increases linearly 
with the light curve's luminosity correction. A, and is well described by cr = 0.055A mag. 
So for a supernova whose luminosity correction A = 0.30 mag, our external distance error 
would amount to less than 0.02 mag. This error is negligible for the majority of observed 
supernovae whose typical |A| < 0.50 mag results in cr < 0.03 mag. The external source of 
error will decrease as the square root of the training set size, which we will be able to expand 
in the future. 
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The time of maximum for each supernova is a non-hnear parameter which cannot be 
solved for analytically in this scheme. It requires an outer iteration to minimize in 
equation (13). For our well-observed light curves, the uncertainty in the time of maximum, 
typically well under one day, has a negligible effect on the parameter errors, with the median 
error increasing by only 10% if the time of maximum is varied ± 1 day. For poorly observed 
light curves whose observations begin ~ 10 days after maximum, the uncertainty in the 
time of maximum increases substantially as does its effect on the parameter errors. We 
discard all light curves beginning more than 10 days after maximum to avoid SN la with 
large uncertainties while maintaining a useful number of objects. When the set of usable 
SN la has grown sufficiently, the precision of this distance indicator might be improved by 
imposing an even stricter requirement for the time of the first observation. 

4 Constructing the Correlation Matrix 

The correlation matrix, C, used in equations (13) and (14) to determine the best param- 
eters and of the light curve fit, is the sum of two parts, C=S-|-N. The noise correlation 
matrix, N, is the correlation matrix for the measurement errors. The signal correlation 
matrix, S, is the correlation matrix which, in the absence of measurement error, estimates 
the expected deviations of the light curves from our model. It is the S matrix which al- 
lows us to use our model despite its known shortcomings. The N matrix is supplied by the 
conscientious observer. 

The S matrix has two parts, diagonal and off-diagonal entries, which, in principle, 
are estimated from our training set. The diagonal entries are estimates of the expected 
deviation of the light curve from the model. We determine the entries along the diagonal 
of each photometric band's block of S by measuring the dispersion of each training set 
member's Ry and Rcoior around the ensemble average Ry and Rcoicn- given in equations (11) 
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and (12) minus each light curve's contribution from measurement error. The result gives, as 
a function of time, the expected errors of noise-corrected data around the best fit model. We 
have plotted this dispersion around the best fit model reconstruction in figure 3 using the 
standard templates (A = 0) as an example. The resulting "grey snakes" comprise 1 sigma 
confidence regions which, when plotted over a best model fit, would be expected to contain 
68% of the data points. The square root of the diagonal entries of the Sv,b,r,i matrix are 
provided in table 2. 
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Figure 3: Dispersion of noisc-free light and color curves around the best fit model. The 
"grey snakes" comprise 1 sigma confidence regions which, when plotted over a best fit model 
fit (here chosen as the A = standard templates), are expected to contain 68% of noise- 
free data points. The square of these functions would be the diagonal entries of the signal 
correlation matrix (S), which, added to the noise correlation matrix (N), would give the 
model correlation matrix (C). To compensate for our inability to determine either data 
covariance (off diagonal elements) or higher order correlations in luminosity we rescale each 
of these functions as described in §4. 
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The off-diagonal entries of the S matrix provide estimates of the correlation between 
model residuals. These correlations are likely to be considerable. There are many more ele- 
ments to consider in determining the point-to-point correlations than we had for the previous 
autocorrelations. While initially we had only to estimate the expected model residuals for 
a given band on a given day of the light curve, the off-diagonal entries require much more 
information. We would need to estimate the amount of covariance between observations on 
different days in the same band, in different bands on the same day, and in different bands 
on different days (In addition, we would have to remove the contribution from measurement 
covariances which observers generally do not supply!). Unfortunately, our sparse training 
set is currently inadequate for quantifying these covariances. A simple two-point correlation 
function suggests that these covariances are most important in B, R, and I but does not 
provide enough information to approximate them adequately. This is the same compromise 
we faced in choosing a simple linear model over one with higher order terms. A detailed 
description of how our data deviates from a linear model requires the same information as 
would be required to establish a more detailed model for our data. Either improvement re- 
quires a larger training set. We echo our minimalist approach to our model with a minimalist 
approach to the S matrix. We employ a diagonal S matrix (with the diagonals determined 
as above) with compensation for the possibility that some photometric bands are more ade- 
quately modeled with our linear model than others. We increase the diagonal elements of the 
signal correlation matrix enough to compensate for our inability to estimate its off-diagonal 
terms. A simple rescaling of the entire signal correlation matrix would be less effective 
because some photometric bands of data (matrix blocks) have more model covariance (i.e. 
non-linear behavior) than others. 

We seek to weight each band's data in the correlation matrix by its abihty to predict the 
parameter in common, A. This approach recognizes that in a linear model some bands may 
be better than others at estimating the parameter A and weights each light curve accordingly. 
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To determine each band's weight, we allow their weights in the correlation matrix to vary 
and maximize the log-likelihood function for the determination of A 

1 " 1 

^^-:,ixl-T.ln-^)- (15) 

^ i=l ^Ai 

Maximizing C is the desired way to determine parameters of the correlation matrix where the 
conventional approach of minimizing would necessarily drive the weights and to zero 
(Rybicki and Kleyna, 1994). A simple exercise shows that maximizing C with respect to the 
weights drives the residuals to the smallest values they can have while still maintaining a 
per degree of freedom near unity. By maximizing C, we optimize our data's ability to predict 
the parameter, A, while simultaneously requiring that the estimated error in A is reasonable. 
Given more training set data, we could estimate the entire correlation matrix by maximizing 
jC. With our current limited training set, we will use our previous prescription to parametrize 
the diagonal blocks of the correlation matrix that correspond to each photometric filter and 
use C to determine the relative weights of those blocks. 

Here C takes the form 

^ = -3(2: —2 -T2 E H-2 z^)) (16) 

^ i=l (^Amlcs + ^Ai„d i=i '^Amlcs + ^A,:„,d 

where the subscripts ind and MLCS denote the value and uncertainty of A as determined 
by independent methods and MLCS respectively. Optimal weighting of the B,V,R, and I 
data minimizes the difference between the independently determined A's (in table 1) and 
the MLCS predicted A's which are a function of the weights in the C matrix. Using equa- 
tion (16) and a downhill simplex method (Press et al 1992) to find its maximum, we have 
determined the relative weights of each band. The diagonal Cv,b,r,i matrix needed for the 
MLCS fit comes from adding the square of the model errors in table 2 to the square of 
the observers' photometry errors, then multiplying the four diagonal blocks of Cv,b,r,i by 
0.37, 11.45,8.85,5.52. These values, determined by maximizing equation (16), indicate that 

28 



the V band hght curves are significantly better predictors of A than B,R, or I data within 
the framework of our linear model. Specifically, a V band observation contributes 5.6, 4.9, 
or 3.9 times as much as a B, R, or I band observation towards determining the luminosity 
correction for a SN la. This result is not astonishing since it is well established that hght 
curves in B,R, and I with differing luminosity can cross each other (Suntzeff, 1993). A cross- 
ing point of such light curves implies the same light or color curve shape for SN la with 
different A's. Such behavior found in B, R, and I diminishes their predictive power in our 
linear model, but the independence of their Ay measurements provide useful estimates of 
extinction. 

We have assumed the correlation matrix for V,B,R, and I data is diagonal, but the 
correlation matrix for V,B-V,V-R, and V-I is certainly not diagonal since an observation of 
my and niB-v made at the same time are anti-correlated with covariance -ay. Formation 
of the correlation matrix for V,B-V,V-R, and V-I data, i.e. Cv,b-v,v-r,v-i, may be done 
directly with care to include the covariance terms of or by means of a simple rotation. 
This involves writing down the rotation matrix. A, which satisfies 



B-V 
V-R 
\V-IJ 



B 
R 
\IJ 



(17) 



and we readily identify 

^V,B-V,V-R,R-I ~ ^ ^V,B,R,I-^- (18) 

Deriving the desired correlation matrix from equation (18) has the advantage of requiring 
the much simpler A and diagonal Cy^g.R,! matrices. 

Meaningful model parameter errors can only come from models which fit the data within 
statistical expectations. We require that both the V light curve model and the color curve 
models give a reduced (per degree of freedom) of 1. To define the confidence region for 
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the distance parameters, we use the covariance matrix of the fit, (L-^C"^L)~^, multiphed 
by the reduced of the hght or color curve fits which measures the parameter of interest. 
This, in effect, is renormahzation of the C matrix, now done on a supernova by supernova 
basis. This renormahzation is not generally a large factor and would presumably become 
unnecesary with a better non-diagonal model for C. 

In the case of the distance modulus error, derived from the visual band data, the error 
is the (1,1) entry of the covariance matrix multiplied by the reduced of the V light curve 
fit: 

al^{L^C-'L)-,\xUV). (19) 

Similarly, the extinction error, as derived from the B-V, V-R, V-I, is the (2,2) entry of the 
covariance matrix multiplied by the reduced of the color curves' fit: 

= {L^C-'L)-,\xUB -V,V-R,V-I). (20) 

The extinction- corrected distance is given by — Ay and its variance is the sum of equations 
(19) and (20) minus twice the covariance of the estimates of ji and Ay: 

^l-Ay = ^2 + - 2{h^C-'h)-,\^ xl{V)xl{B -V,V-R,V-I). (21) 

The extinction-corrected distance error of equation (21) is the previously mentioned 
fitting error. For a particular SN la, its size depends on light curve sampling, measurement 
errors and light curve shape (see §6). These errors provide useful individual estimates of 
distance uncertainty. 

5 Formalized Truncation of Ay 

What is the best way to estimate the absorption by dust given a measurement of excess 
color? We have well-founded a priori knowledge that dust scatters or absorbs light but does 
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not amplify it so the true value of Ay must lie within the range 



0<Av<{n + Mv- miim) 



(22) 



where is a detection limit. Since we measure Ay by dust's reddening effect we say, a 
priori, that dust cannot "blue-en" or brighten an SN la. We could use this knowledge to 
improve our estimate of Ay by truncating any measurement of Ay found to be less than 
zero. Simple truncation carries the disadvantage of improperly treating our useful estimate 
of a Ay from equation (20). What is the best way to use our estimate of Ay and its error 
together with prior knowledge that Ay cannot be less than zero? A straightforward Bayesian 
calculation provides the solution. 

Suppose we were to make an estimate of Ay with value a and normal error aa. Further, 
we have some knowledge of the distribution of the observed Ay, p{Ay). Using Bayes's 
theorem 



p{Av\a,aa) 



p{d\Ay, (7a)p(^y) 
p{a) 



Io^p{Ay)e 



e 




p{Av) 



dAy 



(23) 
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Figure 4: Bayesian estimation of visual extinction by dust. By combining our MLCS mea- 
surement of visual extinction, Ay, and its error with the a priori knowledge that Ay is 
positive we can estimate the likelihood and distribution of the true value for Ay- The top 
panel combines a negative estimate for Ay (=—0.25 ± 0.20) with a gaussian wing (a = 1) 
a priori distribution to yield an estimate for Av{—0.0 ± 0.08). The bottom panel combines 
a positive estimate for Ay (=0.25 ± 0.20) with the same a priori distribution to give an 
estimate for Ay (=0.24 ± 0.16). 
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This "Bayesian filter" provides a probability distribution for the true Ay from which we 
can obtain a best estimate of Ay and its error. A minimalist approach to p{Av) would be 
to assume it is constant over the range in equation (22). This formalized truncation is too 
conservative and unrealistic and we can do better. Since supernovae with very large values 
of Ay are less likely to be part of a sample of detected supernovae, we suggest a reasonable 
form for the observed p{Ay) is a one-sided gaussian which has a maximum at Ay=Q and 
declines for large Ay (see figure 4). We have chosen a—l magnitude of Ay for our p[Ay), 
but we show in §7 that our results are insensitive to the particular value used for a. 

6 Comparing Distances 

The above completes our development of an algorithm to measure extinction-corrected 
distances with multicolor light curve shapes. There are, however, a few complications to 
consider in the practice of measuring distances to supernovae with light curve shapes. 

The "K correction" (Oke and Sandage 1968, Humason, Mayall, and Sandage 1956) 
corrects for the effects of redshift on the measured flux through a fllter of flxed spectral 
response. These corrections can be approximated by measuring the effect of redshifting the 
spectra of SN la supernovae for different redshifts and phases. We can include K corrections 
in our templates given the redshift, or include them in our measurements given the redshift, 
using an assumed time of maximum in the outer iteration. For B and V band data we have 
used the K corrections of Hamuy et al (1993b) and for R and I we have calculated our own 
(see Filippenko et al 1996). The light curves are also affected by time dilation so we contract 
the light curve by (1+z) to return it to the rest frame (Leibundgut et al. 1996, Goldhaber 
et al 1996). 

Using equation (14) we now measure the distance related parameters for a set of well- 
observed SN la assuming that they share the same behavior as our training set. By applying 
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the method as developed for the training set to this independent sample, we compute dis- 
tances for each supernova, and construct a Hubble diagram. Analysis of this diagram shows 
that the MLCS approach gives better precision than the standard candle method. We will 
also compare the results with a sample selected by the same criteria used by Tammann and 
Sandage (1996) to see whether MLCS improves the utility of a "normal" set of SN la. We 
restrict our attention to light curves obtained on a modern photometric system where the 
light curve begins within ten days of maximum light (as determined by our fit). Tests on 
the training set have shown that in order for a light curve to contain luminosity information 
in its shape the first observation must be within ten days of maximum. Every one of the SN 
la in our samples was recorded on digital images and has accurate subtraction of the host 
galaxy background to insure the light curve shapes are free from systematic errors (Boisseau 
and Wheeler 1991). 

Our independent set of twenty supernovae contains ten objects from the Calan/Tololo 
survey (Hamuy et al. 1993a, 1994, 1995, 1996; Maza et al. 1994), two from the hterature 
(Ford ct al. 1993) and eight from our own work (Riess et al. 1996). Table 3 contains SN 
la MLCS parameters. Host galaxy redshifts (column 2) are in the CMB rest frame. The 
hehocentric redshifts (Riess et al 1996, Ford et al 1993, Hamuy 1995) are transformed to 
the Local Group rest frame by the addition of (-30,297,-27) km s~^ in Galactic Cartesian 
coordinates (de Vaucouleurs et al 1991, Lynden-Bell & Lahav 1988). These Local Group 
redshifts are transformed to the CMB rest frame with the addition of (10,-542,300) km s~^ 
(Smoot et al. 1992). For SN 1993ae we have used the redshift of Abell 194 (Chapman 1988) 
of which the host galaxy is a member. For the three SN la with cz < 3000 km s~^ (SN 
1991M, SN 1992G, and SN 1995D) we have corrected their redshifts for their likely infall 
towards Virgo (Schmidt, Kirshner, & Eastman 1992). 

The Galactic extinction measures (column 7) which we use for comparison to our own 
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Table 3: SN la Parameters 
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4.350 


37.79 


0.09 
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0.00 


0.04 


1992P 
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35.50 


0.08 


-0.20 


0.11 


0.00 


1992J 


4.140 


36.75 


0.14 


0.35 


0.18 


0.15 


1991U 


3.991 


35.70 


0.17 


-0.41 


0.75 


0.20 


1991ag 


3.613 


34.15 


0.08 


-0.27 


0.01 


0.13 


1990af 


4.178 


36.84 


0.06 


0.27 


0.00 


0.05 


1992G 


3.299 


32.22 


0.11 


0.02 


0.45 


0.04 


1991m 


3.389 


32.96 


0.10 


0.73 


0.00 


0.01 


1993ac 


3.709 


34.52 
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0.38 


0.00 


0.10 


1994M 


3.859 


35.32 


0.09 


0.21 
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0.00 


1994S 


3.685 


34.24 


0.05 
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1994T 
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4.170 


36.93 


0.22 


0.34 


0.20 


0.45 


1995D 


3.398 


32.76 


0.06 


-0.28 


0.21 


0.00 


1995E 


3.517 


33.76 


0.06 


-0.17 


1.86 


0.00 
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are from Burstein & Heiles (1982). The values of // — Ay and afj,-Av in columns 3 and 4 
are the MLCS extinction-corrected distances and the fitting errors with the former placed 
onto the Cepheid variable distance scale as described in §7 and table 5. The line-of-sight 
extinction estimate, Ay, is listed in column 6 and adding it to the values oi — Ay gives 
the distance estimate without any correction for absorption. Adding these values of /j, to the 
luminosity correction, A, in column 5 gives a distance estimate without correction for either 
the luminosity-light curve relation or absorption. 

In figure 5 wc show the multicolor light curve shape reconstruction for three SN la 
(Riess et al 1996) from our independent sample spanning the range of data quality and 
distance error. For SN 1993ac we have 23 observations beginning shortly after maximum 
light resulting in an extinction-corrected distance error of 0.20 mag. The light and color 
curves of SN 1994Q contain 42 observations beginning shortly after maximum and give an 
extinction-corrected distance error of 0.13 mag. SN 1995D has one of the best sampled 
light and color curves with 107 observations beginning before maximum light and yielding 
an extinction-corrected distance error of 0.06 mag. In general, the size of our predicted 
extinction-corrected distance error depends on the number of observations, the noise in the 
observations, and whether the SN la was first observed before or after maximum light. Half 
of our independent sample of twenty SN la were observed at or before maximum light. 
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SN 1993ac 




Figure 5: Three SN la MLCS fits spanning the range of data quahty in the independent 
sample of twenty SN la. SN 1993ac has 23 noisy observations beginning shortly after maxi- 
mum fight resulting in a distance error of 0.20 mag. SN 1994Q has 42 observations with low 
noise beginning after maximum and has a distance error of 0.13 mag. SN 1995D has 107 
observations with little noise beginning before maximum resulting in a distance error of 0.06 
mag. Photometry is from Riess et al (1996). 
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First, we will use the Hubble diagram as an analytical tool, without reference to a 
distance scale calibration, to determine the precision of our method. Figure 6 shows, for 
comparison, two Hubble diagrams for this independent set of supernovae. In figure 6a we 
have fit the best light curve shape to each supernova to estimate the distance without any 
correction for intrinsic luminosity variation or extinction. In figure 6b we have plotted the 
MLCS extinction-corrected distances to each SN la which accounts for intrinsic luminosity 
variation and extinction (see table 4) . The reduction in dispersion is dramatic. The improve- 
ment in distance precision comes from deriving the correlation between luminosity, and light 
and color curve shape from our training set of SN la, and then applying these relations to the 
independent sample. Because the training set and the independent set have no overlapping 
members, the reduction in dispersion is a powerful demonstration of the effectiveness of the 
MLCS method. 

Table 4 compares the distance estimates for different assumptions by measuring the 
dispersion and on the Hubble diagram. In each case we have made a custom reconstruction 
of the light and color curves. The first row is the "standard candle" assumption for which we 
disregard the light curve shape-luminosity correction. A, by adding it back to the distance, 
// -|- A, and make no correction for extinction. Next, we use the MLCS distance modulus, fj, 
which includes the luminosity information, but make no allowance for absorption. Following 
this we include a correction for only the Galactic component of extinction by using Burstein 
& Heiles (1982) absorption measures. Finally, we use the full MLCS method to estimate 
the extinction-corrected distance. We compare the different methods for three different 
subsamplcs: all twenty SN la, all SN la minus the single most highly reddened object (SN 
1995E), and with a color cut (Vaughan et al 1995, Tamman & Sandage 1995, Hamuy et al 
1995) which discards objects outside the color range —0.25 < {B — V)max < 0.25, leaving 
eighteen objects. 
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DISTANCE MODULUS 

Figure 6: Hubble Diagrams for SN la with velocities in the COBE rest frame on the Cepheid 
distance scale (Sandage et al 1994, 1996). All velocity errors are 300 km reflecting a plau- 
sible estimate of random velocities with respect to the Hubble flow, (a) Distances estimated 
with a standard luminosity assumption and no correction for extinction. This method yields 
(7^=0.52 and Ho=52 ± 8 (statistical) km s~^ Mpc~^(b) Distances from the MLCS method 
which makes a correction for intrinsic luminosity variation and total extinction as deter- 
mined from the light and color curve shapes. This method yields ay—0.12 and Ho=65 ± 3 
(statistical) km s~^ Mpc~^. 
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Table 4: SN la distance comparisons (CMB frame, a-^ei — 300 km s ^) 



Correction 


Form 


all SN la, N=20 


w/o 


SN1995E, N=19 


color cut, N=18 


luminosity 


Av 




a Pr 


a 


Pr 


a Pr 


none 


none 


H+A 


0.52 


0.40 




0.33 


MLCS 


none 




0.46 0.015 


0.20 


< 1x10-6 


0.20 2x10-5 


MLCS 


galactic 


— Av,gal 


0.46 0.25 


0.17 


0.018 


0.16 0.0085 


MLCS 


MLCS 




0.12 1x10-^ 


0.13 


7x10-'' 


0.13 0.00013 



Moving down table 4, each successive row represents a refinement in our distance mea- 
suring technique. The result of successive improvements in the method can be seen in the 
decreasing dispersion on the Hubble diagram. Including the heavily reddened SN 1995E in 
the sample demonstrates the power of MLCS in dealing correctly with reddened objects but 
masks the gradual improvement in distance to be made from various aspects of the method. 
It is important to note that even the distances of the color cut sample can be improved with 
MLCS. This shows that the distance precision of "normal" or "Branch-normal" SN la can 
be significantly enhanced using hght and color curve information. The MLCS method is not 
just "reining in" the extreme SN la, but rather is improving the distance measures to all of 
the SN la. 

To demonstrate the statistical significance of each level of improvement, we have in- 
cluded in table 4 the probabilities, Pr, that the observed improvement in dispersion could 
occur from a random set of distance corrections. These probabilities give the likelihood for 
the null hypothesis; that our distance "corrections" have no relation to the true SN la dis- 
tances. In a Monte Carlo simulation we apply a set of random corrections chosen from the 
distribution of proposed corrections and see how often the dispersion is as low or lower than 
the actual improved dispersion. The value of Pr gives the probability that the observed (or 
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a greater) decrease in dispersion for each successive distance refinement occured by chance. 
The results show MLCS luminosity and extinction corrections are highly significant regard- 
less of the sample criterion. Both the MLCS luminosity correction and extinction correction 
strongly reject the null hypothesis that they are unrelated to the true distance of the SN la. 
Even the small corrections for galactic extinction do more good than harm. Yet, accounting 
only for the Milky Way's contribution to the total absorption fails to account for host galaxy 
extinction which in a few cases can be substantial. Using all the predictive power of MLCS 
gives remarkably low values for the dispersion with an exceedingly small probability that this 
improvement in dispersion occured by chance. With a conservative estimate for the peculiar 
velocity associated with each field galaxy of 300 km s~^ our observed dispersion of 0.12 mag- 
nitudes implies a typical distance precision of 5%. The improvement in distance precision 
by including a correction for host galaxy extinction with the conventional reddening law is 
the first demonstration that such corrections can be successfully made. 

Our Monte Carlo simulation demonstrates that even for a set of SN la selected by 
color, such as those used by Sandage et al (1996), the MLCS method makes a significant 
improvement in the precision of the distances. For the color cut sample, the probability that 
both our luminosity and extinction corrections would improve the dispersion from 0.33 mag 
to 0.13 mag by chance is less than one in a million. 

Both SN 1992K and SN 1995E provide instructive examples of how this method leads 
to improved distance estimates. Both objects appear to be dim, which for a standard candle, 
suggests they are at a great distance. Assuming a standard candle luminosity places each 
of them much further away than their redshift implies (see figure 6a). Yet, there are clues 
in the light and color curves of these objects which indicate that they are dim for different 
reasons. The rapidly declining V light curve of SN 1992K and the color evolution of its 
B-V curve are nearly identical to the photometric behavior of the subluminous SN1991bg, a 
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member of the training set (Hamuy et al 1994). Application of the MLCS method estimates 
SN 1992K to be A =1.25 magnitudes dimmer than the standard SN la, though its Ay is 
only 0.01 mag. This correction to the luminosity is independent of the SN's redshift since it 
depends only on the light curve shape, but it is reassuring to note that accounting for this 
object's intrinsic faintness shifts its distance quite precisely onto the Hubble line. 

For the case of SN 1995E, the dim appearance that places this SN la below and to the 
right of the Hubble line is not intrinsic to the supernova, but rather, is a result of absorption, 
as shown by MLCS. This supernova was found on the spiral arm of NGC 2441. The shape of 
its light and color curves suggests a fair resemblance to the standard SN la event (A = —0.17 
mag), but all its measured colors are systematically displaced to the red, as would occur from 
absorption by dust. By fitting the shape of the color curves independently from the value of 
the color, we can measure the color excesses. Assuming a standard reddening law (see Riess, 
Press, and Kirshner 1996b) we estimate the visual band extinction to be 1.86 magnitudes. As 
in the case of SN 1992K, correcting the luminosity of SN 1995E makes its distance consistent 
with its redshift measurement. Using a color cut requiring —0.25 < {B — V)max < 0.25, both 
of these objects would be discarded. Yet, we can keep these objects (and others like them) 
in the sample and use the MLCS method to distinguish between supernovae which are 
intrinsically dim and those which are dimmed by dust absorption. 

We can make two significant checks of our extinction-corrected distances by examining 
the fit of the Hubble line to the data. We examine the linearity of the Hubble law by 
measuring the slope of the relation between the MLCS distance modulus and the logarithm 
of the redshift. Assuming that space is Euchdean for our modest redshifts, the expectation 
of this slope is 0.2. Lauer & Postman (1992) found a slope of 0.1992 ± 0.006 using brightest 
cluster galaxies and Jcrjen & Tammann (1993) found a slope of 0.1988 ± 0.006 using the 
mean of a number of distance indicators to fifteen clusters. Using the extinction-corrected 
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distances and errors of table 3 yields a slope of 0.2010 ± 0.0035 which is consistent with 0.2 
and with the two previous results. The small error on our slope over the distance interval 32.2 
< < 38.3 makes this the most precise check of this classical test of cosmology. Finally, 
we can examine the goodness of fit of our Hubble hne to the extinction-corrected MLCS 
distances. The of the fit using the independently determined distance errors in table 3 
and a^ei = 300 km s~^ is 13 for 19 degrees of freedom, which is within the expectation of x^- 
The value of is strongly dependent on the assumed random velocity of the field galaxies 
hosting our SN la due to our low dispersion and small distance errors. A random velocity 
error for our field galaxies in the range 125 km s~^ < a^ei < 300 km s~^ gives a within its 
likely range of 13 to 25 and is consistent with other determinations of this random velocity 
component (Marzke 1995, Davis & Peebles 1983). All indications suggest that the MLCS 
method provides remarkably precise, extinction-corrected distances which are a significant 
improvement over SN la distances from previous methods. 

7 Discussion 

Our intent has been to describe the MLCS method in enough detail so others can take 
advantage of the precise distance estimates it provides. The strength of MLCS lies in its 
ability to disentangle the effects of absorption and intrinsic luminosity variation while provid- 
ing meaningful error estimates. These measures are derived from the distance independent 
observables of multicolor light curve shapes. The accuracy of the MLCS relative distance 
measures has been well established on an independent set of twenty SN la on the Hubble 
diagram. 

To place our MLCS distances on the established absolute distance scale we use the 
luminosity calibration for a number of SN la with an independent distance indicator of high 
precision. At present, there are three SN la whose hght curves meet our MLCS quality criteria 
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Table 5: SN la distance calibration 



SN la 


IJ'MLCsicf) 


/J'Cephio') 


fJ'MLCS — jJ'Cephify) 


A(a) 


Mv - Mv,81b{(J) 


1972E 
1981B 
1990N 


27.21*(0.09) 
30.38(0.07) 
31.04(0.14) 


28.08*(0.10) 
31.10(0.20) 
32.00(0.23) 


0.87(0.13) 
0.72(0.21) 
1.01(0.0.27) 


-0.33(0.04) 
0.01(0.04) 
-0.27(0.04) 


-0.28(0.21) 
0.00(0.00) 
-0.07(0.23) 



* distance modulus uncorrected for extinction 



(modern photoelectric photometry with observations less than ten days after maximum) 
and whose distances have been measured with Cepheids observed with the Hubble Space 
Telescope. 

The three SN la, SN 1972E, SN 1981B, and SN 1990N, are hsted in table 5: column 
(2) gives the extinction-corrected MLCS distances on the distance scale of SBF-PNLF-TF 
(table 1), column (3) gives best distances as determined by HST Cepheid measurements 
(Sandage et al 1994, 1996), column (4) gives the differences between columns (3) and (2), 
column (5) gives the MLCS luminosity correction, and column (6) gives the Sandage et al 
(1996) My minus the My for the "standard" shaped SN 1981B. By comparing our precise 
MLCS relative distances to the trustworthy Cepheid absolute distances we see a consistent 
difference of 0.81±0.10 mag. Adding this difference to all of our distance estimates places our 
distance indicator on the Cepheid distance scale. Further, it provides an absolute luminosity 
cahbration for the standard (A = 0) SN la of My = —19.36 ± 0.10 mag at B maximum. 
This calibration for the standard My template is used in table 2. The extinction-corrected 
distances listed in table 3 include this offset which places them on the Cepheid distance 
scale. We can calculate the Hubble constant by fitting 

logv = .2(/x - Ay) + log Hq - 5 (24) 

to the distances and velocities in table 3 using the tabulated errors a^^Av assumed 
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Table 6: SN la Hubble constant comparisons (CMB frame, a-vei — 300 km s ^) 



Correction 


Form 


all SN la, N=20 


w/o SN1995E, N=19 


color cut, N=18 


luminosity 


Av 




Ho 


Ho 


Ho 


none 


none 


H+A 


52 ±8 


54 ±6 


55 ±5 


MLCS 


none 




61 ±9 


63 ±4 


62 ±4 


MLCS 


galactic 


1^ — Av,gal 


63 ±9 


65 ±4 


64 ±4 


MLCS 


MLCS 




65 ± 3 


65 ± 3 


65 ± 3 



random velocity of 300 km s~^. The result is a Hubble constant of 66 ± 3 km s~^ Mpc~^ 
where the uncertainty is internal and incorporates a 0.03 mag uncertainty in the Hubble 
line and a 0.10 mag uncertainty in the placement onto the Cepheid distance scale. We also 
see that the MLCS luminosity corrections (column 5) are consistent with the luminosity 
differences as determined by the Cepheids (column 6). 

For comparison we perform the same calculation for the other distance and extinction 
methods using the three SN la samples. We apply each method to the cahbration set to 
determine the offset onto the Cepheid scale, and to the distant set to determine the Hubble 
line. For the methods which do not provide their own error estimates, we assume a constant 
distance error of the size necessary to get the expected chi'^ of the fit to the Hubble line. In 
table 6 we list the determinations of the Hubble constant and internal error for the different 
methods and samples. 

The color cut sample provides us with a set of SN la which are directly comparable 
to the "normal" set used by Sandage et al (1996). Using the assumption of a standard 
luminosity and the color cut sample gives i^o = 55 ± 5 km s~^ Mpc~^ which is consistent 
with the Sandage et al (1996) values of Ho{B) = 56 ± 4 and Ho{V) = 58 ± 4 km s"^ Mpc^^ 
obtained with the same standard brightness assumption and sample criteria, though with a 
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different set of light curves. 

In §6 and table 4 we demonstrated that each of the distance measuring improvements 
en route to the complete MLCS method is highly significant and should be employed to 
get the best result. Using the MLCS method to measure extinction-corrected distances and 
errors yields a Hubble constant of i^o = 65 ± 3 km s~^ Mpc~^ or i^o = 63 ± 3 km s~^ Mpc~^ 
using constant weighting. The only significant change in the Hubble constant arises from 
including the luminosity correction, A. This is because the mean luminosity of the three 
nearby calibrators is ~ 15% dimmer than the mean luminosity of the twenty SN la in the 
distant sample. That these nearby SN la are dimmer than the distant SN la does not imply 
an anti-selection bias because our sample is neither volume nor magnitude limited. If our 
sample was complete in volume or magnitude, we would observe a tremendous increase in 
the number of observed SN la at large distances. Figure 6 shows this is clearly not the case. 
We elaborate on this point below. 

A complete error budget for any of the values of Hq in table 6 would consist of the 
stated internal error added in quadrature to the estimated uncertainty in the Cepheid zero 
point, roughly 0.15 mag (Feast and Walker). For the MLCS extinction-corrected distances, 
this gives a value of ifo = 65 ± 6 km s~^ Mpc~^. 

The MLCS distances to SN la in the Hubble flow combined with redshifts provide 
the necessary information to estimate the peculiar velocity component of the each objects' 
radial motion. Plotted on the sky, the velocity residuals show a dipole pattern indicative 
of the motion of the Local Group with respect to a frame deflned by the supernovae. Our 
preliminary analysis of this motion with a subset of 13 SN from the current set and without 
our MLCS reddening information was consistent with convergence to the cosmic microwave 
background frame and inconsistent with the Lauer & Postman frame (1994) at 7000 km 
s~^ (Riess, Press, Kirshner 1995b). The same analysis performed at higher precision using 
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the independent set of twenty SN la and the MLCS method (which now includes extinction 
corrections) yields an even stronger detection of the Local Group motion with similar results. 
In the future, when the sample of SN la has grown, we will revisit our analysis of the Local 
Group motion from SN la. 

Sections 2-5 outline the MLCS technique in sufficient detail to allow for future redeter- 
mination of the necessary templates and functions as more supernovae become available or 
are deemed desirable to include in the training set. The training set for this paper uses all 
supernovae for which accurate and extensive photoelectric photometry is currently available 
as well as definitive relative distance estimates and Ay estimates. These included from the 
Phillips (1993) set, SN 1980N, SN 1981B, 1986G, 1989B, 1990N, 1991T, 1991bg, and 1992A. 
We exclude SN 19711 for which only photographic photometry is available . We also use SN 
1994ae (Riess et al 1996) for which there are detailed photometric light curves and an inde- 
pendent distance (Dell' Antonio 1995). Ideally we would use only SN la's with no evidence 
of extinction, but currently we lack a sufficient number of training set objects to discard any. 
Instead, we estimate Ay from the color differences of unreddened SN la's of similar light 
curve shape and propagate the resulting uncertainty in Ay to our distance errors. 

Richmond et al (1995) have discussed a conspicious inconsistency with a preliminary 
SBF distance to NGC 4526 (John Tonry, private communication) and the distance from 
the shape of the SN 1994D hght curve. The difference of 0.7 magnitudes for the SBF and 
supernovae distance moduli to NGC 4526 imply a "greater than 3-sigma" mutual rejection. 
The case can be simplified by comparing the results for SN 1992A in NGC 1380 in Fornax to 
SN 1994D. These supernovae have virtually identically shaped hght curves, but SN 1994D 
peaked 0.65 magnitudes brighter than SN 1992A implying that it was 33% closer. SBF 
suggests these two supernovae arc in galaxies at similar distances (Phillips 1993). Current 
HST Key project observations of Cepheids in Fornax may eventually help to "break the tie" . 
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Inclusion of SN 1994D has the effect of decreasing the precision of the distance estimates 
for the independent sample of objects, but the training set is large enough that the effect 
is not substantial. The dispersion of distances on the Hubble diagram for MLCS increased 
from 0.12 mag to 0.22 mag and the increased from 13 to 39. If SN 1994D were truly at 
the SBF distance and was representative of objects in the independent set we would expect 
to see no significant change in the Hubble diagram for including this object. If however, 
the SBF distance is incorrect or this one object is a special case, then we would expect to 
see a substantial change in the value of as we have. A final possibility would be that 
SN 1994D is represent it ive of a departure from the behavior of other SN la and it has only 
been due to chance that we have not seen any others in the independent sample. In this 
case, the increased dispersion would be a more accurate measure of the true dispersion of 
well-behaved SN la's and "94D-like" objects treated as a whole. Future discovery of such 
"SN 1994D-like" objects could lead to a new class of objects at least for the purposes of 
distance measurement. Until either another such object is observed or the SBF distance is 
finalized or an independent distance is measured we remain agnostic. 

In section 5 we developed a formal way to combine our measurement of Ay with our a 
priori understanding of dust. This method requires some description of the distribution of 
Ay values for supernovae discovered in galaxies. The most conservative estimate for p{Av) is 
that it is constant over the range of equation (22). This amounts to truncation of extinctions 
less than zero. A gaussian wing which is a maximum at Ay— 0.0 and is parameterized by 
its second moment seems more plausible. Either assuming a constant value for p{Av) or 
a gaussian wing with 0.5 < a < oo yields a high level of distance precision. The Hubble 
diagram dispersion is insensitive to the value of a over this range, and the Hubble constant is 
insensitive to any particular parameterization of p{Av) including requiring Ay—O (see table 
6). Using a p{Av) with a < 0.5 is too restrictive and amounts to discarding any extinction 
information. Determination of the best value for a by maximizing the log-likelihood of 
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equation (15) on the Hubble diagram yields a— 0.9 mag. It is important to consider sample 
selection effects when choosing a form for the observed p{Ay). Clearly the width of p{Av) 
could be a function of the search characteristics since these determine how likely it is to 
find an SN la obscured by dust. Monte Carlo calculations provide an estimate of the effect 
of misjudging p{Av), and with enough data, one could easily solve for the distribution of 
observed Ay- 

Tammann and Sandage (1995) have questioned the validity of our training set distances; 
"Supernovac used by Phillips (1993) ...based on distances determined by the Tully-Fisher 
and the surface brightness fluctuation methods... clearly deviate from the [Cepheid calibrated] 
supernovae values." This criticism is misdirected. While it is true that the absolute distance 
scale as determined by Cepheid variables is not necessarily consistent with the Tully-Fisher 
and SBF distances, we only employ relative distances of these methods to measure the relative 
luminosity variation of SN la out to the Virgo cluster (Kennicutt, Freedman, & Mould 1995). 
The selection bias pointed out by Sandage (1988a,b 1994a,b) to distort Tully-Fisher distances 
does not apply to the distances in our training set which were derived from the bias-corrected 
inverse Tully-Fisher relation (Strauss & Willick 1995, Schechter 1980, Pierce 1996). 

Tammann and Sandage (1995) have also commented on the statistical problem raised 
by a light curve-luminosity relation stating that "the Malmquist effect on the distant SN la 
would be overwhelming". Specifically, Tammann and Sandage (1995) argue that the distant 
SN la occupy "a volume of about 3 x 10^ larger, on average, than the volume of the three 
nearby local calibrators" . Therefore, one would expect a substantial selection bias favoring 
brighter SN la in the distant sample as compared to the SN la nearby. In fact, the mean 
SN la luminosity for the distant sample we use is 0.34 magnitudes dimmer than for the set 
of three nearby calibrators. Is this surprising? Not very. The selection bias, as stated by 
Tammann and Sandage (1995), is estimated for a complete sample of SN la with a well- 
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defined limiting magnitude. This is clearly not case for our sample or for theirs. Inspection 
of the Hubble diagrams of figure 6 (or of Tammann & Sandage's) does not show the number 
of SN la increasing with distance as 10^-^^ as expected for a complete search of increasing 
volume. This point can be made quantitatively by performing a simple (y^) test (Schmidt 
1968); which has an expectation value of 0.5 for a uniformly distributed sample. For the 
distant set of 20 SN la { y^ )=0.09, which shows that this sample is concentrated nearby. 
Even limiting the test to the set of ten SN la discovered during the uniform Calan/Tololo 
survey yields (^^-^)=0.16. Real samples of observed SN la are not distributed the way 
assumed by Tamman & Sandage (1995). The chances of the three local calibrators being 
0.34 magnitudes dimmer than the independent sample is 6 %. This is about as likely as 
randomly picking the winner of the next baseball American League Championship series 
(ALCS). 

An interesting question to consider is what explosion or progenitor parameters could 
explain our empirically determined variation in light curve shape, luminosity and color? 
Current models have attempted to explain the inhomogenity of supernovae in one of two 
ways. Hofiich, Khokhlov, & Wheeler (1995) and Hoflich & Khokhlov (1996) have found 
that a variation in the density at which the defiagration burning front transitions to a 
detonation wave affects the amount of ^^Ni produced in the explosion. A late transition 
gives the outer layers time to preexpand resulting in a small nickel production. The reduced 
nickel heating diminishes the temperatures in the expanding envelope and photosphere. This 
results in rapidly dropping opacities. Consequently, the photosphere recedes fast and the 
stored energy is emitted over a short period of time. Conversely, an early turnover of the 
deflagration into a detonation front results in a large amount of nickel which produces a bright 
and hot supernova whose opaque layers keep the radiative energy loss comparably low. This 
mechanism works for deflagrations, delayed detonations and pulsating delayed detonations. 
Qualitatively, the observed correlation between luminosity and the photometric parameters 
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are reproduced (Hoflich & Khokhlov 1996). 

Another theoretical approach to matching the observations involves exploding progen- 
itors of varying and generally sub-Chandrasekhar mass. A layer of helium accumulates at 
the base of the static hydrogen-burning zone. Sudden burning of this helium at the base of 
the layer sends a shock that can trigger a carbon detonation at the interface. If this does not 
occur, a second chance at carbon detonation can come when the shock propagates around 
the star and converges on the opposite side. The variation in progenitor mass suggests a 
simpler connection to the amount of nickel produced than the Chandrasekhar models. Again 
the variation in nickel yield is expected to match the observed correlations in supernova ob- 
servables. These models have been successful in one and two dimensions, but it remains to 
be seen if the abundance and the velocity distribution of the intermediate mass elements 
produced from nucleosynthesis matches the observations (Livne and Arnett 1995, Woosley 
& Weaver 1994, Livne 1990, Livne & Glasner 1990) 

Finally, improvements which can be made in the MLCS technique when we acquire a 
larger training set of supernovae. First, the "grey snakes" of figure 3 and the relative weights 
of the B,V,R, and I data determined from equation (16) suggest that greater precision could 
be attained by including a quadratic term in the light curve shapes model for B,R, and I 
band data. An additional order in the fit would demand a larger training set of data to avoid 
overfitting the details of the training set objects. A larger training set would also support 
the determination of a more detailed correlation matrix with estimates of model residual 
covariances. 

The MLCS technique provides an exceedingly precise way of measuring extinction- 
corrected distances with uncertainty estimates. The redshifts to which Type la supernova 
light curves can measure distances have interesting implications for cosmological measure- 
ments. Two teams are currently laboring to find and measure SN la at 0.3 < z < 0.6 
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(Perlmutter et al 1996, Schmidt et al 1996) with the intent of measuring the cosmological 
deceleration of the Universe, go- At these redshifts SN la hght curves are difficult to obtain, 
and informative spectra even harder, so it is sensible to use all the available SN la. The 
uncertainty in go is proportional to the distance precision divided by the square root of the 
number of objects. The MLCS method should help with both of these factors by including all 
well-observed SN la in the sample, and increasing the precision of SN la distance measures. 

We are again grateful to Mario Hamuy, Mark Phillips, Nick Suntzeff and the entire 
Calan/Tololo collaboration for the opportunity to study their exceptional data before pub- 
lication. We have greatly benefited from discussions with Brian Schmidt, George Rybicki, 
and Peter Hofiich. This work was supported through NSF grants AST 92-18475 and PHY 
95-07695. 
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